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\ Semileptonic and nonleptonic decays of the B c meson to B s and B mesons, caused 

by the c — > s, d quark transitions, are studied in the framework of the relativistic 

00 ' 

o 

spectator b quark is used to simplify calculations while the final s and d quarks 
in the B s and B mesons are treated relativistically. The decay form factors are 
explicitly expressed through the overlap integrals of the meson wave functions in 
the whole accessible kinematical range. The obtained results are compared with the 
predictions of other approaches. 
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I. INTRODUCTION 

The B c meson discovered at Fermilab [1] is the only quark- ant iquark bound system (be) 
composed of heavy quarks (b, c) with different flavours, thus being flavour asymmetric. The 
investigation of the B c meson properties (mass spectrum, decay rates, etc.) is therefore of 
special interest compared to symmetric heavy quarkonium (bb, cc) ones. The difference of 
quark flavours forbids the annihilation of B c into gluons. As a result the pseudoscalar be 
state is much more stable than the heavy quarkonium one and decays only weakly. It serves 
as a final state for the pionic and radiative decays of the excited be states (lying below the 
BD threshold). Experimental study of the B c mesons is planned both at the Tevatron and 
Large Hadron Collider (LHC) (for a recent review see e.g. [2] and references therein). 

Since both quarks (b, c) composing the B c meson are heavy, their weak decays contribute 
comparably to the total decay rate. Thus there are two distinctive decay modes: (i) b c,u 
with c quark being a spectator, and (ii) c — > s, d with b quark being a spectator. The 
transition (i) induce the semileptonic B c decays to charmonium and D mesons, while the 
transitions (ii) cause the B c decays to B s and B mesons. The kinematical ranges of these 
semileptonic decay modes are substantially different. The square of the four momentum 



transfer to the lepton pair extends from to g 2 iax ~ 10 GeV 2 for the decays to charmonium 
and g^ ax ~ 18 GeV 2 for decays to D mesons, but only to g 2 nax ~ 1 GeV 2 for decays to B and 
B s mesons. Thus the kinematical range for the decay mode (i) is appreciably larger than for 
the decay mode (ii). Otherwise in the B c rest frame the maximum recoil three momentum 
of the final charmonium and D meson turns out to be of order of their masses, while that 
of final B and B s mesons is much smaller than the meson masses. 

The weak B c decays to charmonium and D mesons were studied at length in our recent 
paper [3]. Here we consider the weak B c decays to B s and B mesons within the relativistic 
quark model. The model is based on the quasipotential approach in quantum field theory 
and was fruitfully applied for describing the electroweak decays and mass spectra of heavy- 
light mesons, heavy quarkonia [4-9] and B c meson [10]. The relativistic wave functions 
obtained in the latter paper are used below to calculate the transition matrix elements. 
The consistent theoretical description of B c decays requires a reliable determination of the 
q 2 dependence of the decay amplitudes in the whole kinematical range. In most previous 
calculations the corresponding decay form factors were determined only at one kinematical 
point either q 2 = or q 2 = g 2 ^ and then extrapolated to the allowed kinematical range 
using some phenomenological ansatz (mainly (di)pole or Gaussian). Our aim is to explicitly 
determine the q 2 dependence of form factors in the whole kinematical range thus avoiding 
extrapolations and reducing uncertainties. 

The paper is organized as follows. In Sec. II we describe the underlying relativistic 
quark model. The method for calculating matrix elements of the weak current for c — > s, d 
transitions in B c meson decays is presented in Sec. III. Special attention is paid to the 
dependence on the momentum transfer of the decay amplitudes. The B c decay form factors 
are calculated in the whole kinematical range in Sec. IV. The q 2 dependence of the form 
factors is explicitly determined. These form factors are used for the calculation of the B c 
semileptonic decay rates in Sec. V. Section VI contains our predictions for the energetic 
nonleptonic B c decays in the factorization approximation, and a comparison of our results 
with other theoretical calculations is presented. Our conclusions are given in Sec. VII. 
Finally, the Appendix contains complete expressions for the decay form factors. 



In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation [11] of the Schrodinger 



II. RELATIVISTIC QUARK MODEL 
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where the relativistic reduced mass is 
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and Ei, E 2 are the center of mass energies on mass shell given by 
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Here M = E x + E 2 is the meson mass, m 12 are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 
reads 

2 [M 2 - (mi + m 2 ) 2 ][M 2 - (m, - m 2 ) 2 ] 

4M 2 ' 1 J 

The kernel V(p, q; M) in Eq. (1) is the quasipotential operator of the quark- ant iquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [4] 

V(p, q; M) = fZ!(p)iZ 2 (-p)V(p, q; M) Ul (q)u 2 (-q), (5) 

with 

V(p, q; M) = ^ s ZV( k K72 + KL(k)r?r 2;M + V^ nf (k), 
where a s is the QCD coupling constant, D^ u is the gluon propagator in the Coulomb gauge 

Air 4-ir ( PfrA 

£°°(k) = -^, £,«( k ) = \& - ^-J , D« = J D" = > (6) 

and k = p — q; 7 M and u(p) are the Dirac matrices and spinors 
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Here <r and x are the Pauli matrices and spinors; e(p) = ^p 1 + m' 2 . The effective long-range 
vector vertex is given by 

r„(k) = 7^ + ( 8 ) 

where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

V v (r) = (1 - e)Ar + B, 

V s (r) = eAr, (9) 

reproducing 

V COQ{ (r) = V s (r) + V v (r) = Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v 2 /c 2 without 
and with retardation corrections to the confining potential, can be found in Refs. [4] and 
[5, 10], respectively. The structure of the spin-dependent interaction is in agreement with 
the parameterization of Eichten and Feinberg [13]. The quasipotential for the heavy quark 
interaction with a light antiquark without employing the expansion in inverse powers of the 
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FIG. 1: Lowest order vertex function contributing to the current matrix element (11). 



light quark mass is given in Ref. [6]. All the parameters of our model like quark masses, 
parameters of the linear confining potential A and B, mixing coefficient e and anomalous 
chromomagnetic quark moment k are fixed from the analysis of heavy quarkonium masses [4] 
and radiative decays [7]. The quark masses nib = 4.88 GeV, m c = 1.55 GeV, m s = 0.50 GeV, 
m u rf = 0.33 GeV and the parameters of the linear potential A = 0.18 GeV 2 and B = —0.16 
GeV have usual values of quark models. The value of the mixing coefficient of vector and 
scalar confining potentials e — — 1 has been determined from the consideration of the heavy 
quark expansion for the semileptonic B — > D decays [8] and charmonium radiative decays 
[7]. Finally, the universal Pauli interaction constant k — — 1 has been fixed from the analysis 
of the fine splitting of heavy quarkonia 3 Pj- states [4]. Note that the long-range magnetic 
contribution to the potential in our model is proportional to (1 + k) and thus vanishes for 
the chosen value of k — — 1. It has been known for a long time that the correct reproduction 
of the spin-dependent part of the quark- ant iquark interaction requires either assuming the 
scalar confinement or equivalently introducing the Pauli interaction with k — — 1 [4, 5, 14] 
in the vector confinement. 



III. MATRIX ELEMENTS OF THE ELECTROWEAK CURRENT FOR c -»• s, d 

TRANSITIONS 

In order to calculate the exclusive semileptonic decay rate of the B c meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current Jjf = <?7 M (1 — 75)0, 
associated with c — > q (q = s or d) transition, between a B c meson with mass Mb c and 
momentum pb c and a final meson F (F = B^ or £?M) with mass Mp and momentum pp 
takes the form [15] 

(F(p F )\J™\B c (p Bc )) = J ^^^ PF (p)r,(p,q)^ cPflc (q), (11) 

where T^(p,q) is the two-particle vertex function and $m Pm are the meson (M = B C ,F) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum p«. 

The contributions to T come from Figs. 1 and 2. The contribution T < - 2 ' ) is the consequence 
of the projection onto the positive-energy states. Note that the form of the relativistic 
corrections resulting from the vertex function is explicitly dependent on the Lorentz 
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FIG. 2: Vertex function taking the quark interaction into account. Dashed lines correspond to 
the effective potential V in (5). Bold lines denote the negative-energy part of the quark propagator. 



structure of the quark- ant iquark interaction. In the leading order of the v 2 /c 2 expansion 
for B c and in the heavy quark limit mj, — > oo for B S ,B only contributes, while r^ 2 ) 
contributes already at the subleading order. The vertex functions look like 



and 



r«(p, q) = u q [p q ) llx {\ - 7 5 )M c (g c )(2vr) 3 ( 5(p fc - q b ), 



If)(p,q) = uMuM^l - 7i) 6c( ^^ ) ( ^ ) 7?V(p b - q 6 ) 

+V(p b - qfc) e + e ( g j TiW 1 - l 5 i)}u c (q c )u b (q b ), 



(12) 



(13) 



where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, k = p q — A; k' = 
q c + A; A = p F -p Sc ; 



A( _) = e(p) -(m 7 + 7 (7p)) 
1 J ~ 2e(p) 



Here [15] 



P,,6 = ^ b (p)^-±J2n ( - l \p F )p\ 



M B C i= i 



= 1 
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and are three four-vectors given by 



p' 



p l pi 



M(E + M) 



, E = ^p 2 + M 2 . 



It is important to note that the wave functions entering the weak current matrix element 
(11) are not in the rest frame in general. For example, in the B c meson rest frame (p# c = 0), 
the final meson is moving with the recoil momentum A. The wave function of the moving 
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meson \&i?A is connected with the wave function in the rest frame ^fo = by the 
transformation [15] 

*fa( P ) = Dl/\RYjDi /2 (RYj* F0 (p), (14) 

where R w is the Wigner rotation, La is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D l l 2 (R) in spinor representation is given by 

(J °^D 1 q / ;(RYj = S-\ Pqtb )S(A)S(p), (15) 

where 

/ e(p) + m / ap \ 

5(p) = V^^( 1 + ^)T^J 

is the usual Lorentz transformation matrix of the four-spinor. 

The general structure of the current matrix element (11) is rather complicated, because it 
is necessary to integrate both with respect to d 3 p and d 3 q. The 5-function in the expression 
(12) for the vertex function permits to perform one of these integrations. As a result 
the contribution of to the current matrix element has the usual structure of an overlap 
integral of meson wave functions and can be calculated exactly (without employing any 
expansion) in the whole kinematical range, if the wave functions of the initial and final 
mesons are known. The situation with the contribution r^ 2 ^ is different. Here, instead of a 
5-function, we have a complicated structure, containing the potential of the gg-interaction 
in a meson. Thus in the general case we cannot accomplish one of the integrations in the 
contribution of r^ 2 ^ to the matrix element (11). Therefore, one should use some additional 
considerations in order to simplify calculations. The main idea is to expand the vertex 
function T^ 2 \ given by (13), in such a way that it will be possible to use the quasipotential 
equation (1) in order to perform one of the integrations in the current matrix element (11). 

The natural expansion parameters for B c decays to B s , B mesons are the active c and 
spectator b quark masses. However, the heavy c quark undergoes the weak transition to the 
light s or d quark. The constituent s, d quark masses are of the same order of magnitude as 
the relative momentum and binding energy, thus we cannot apply the expansion in inverse 
powers of their masses. The heavy quark expansion in l/m c ^ significantly simplifies the 
structure of the r^ 2 ^ contribution to the decay matrix element, but the momentum p depen- 
dence of the light quark energies e q (p) still prevents to perform one of the integrations. It is 
important to note that the kinematically allowed range for B c decays to B s and B meson is 
not large (| A max | = (M| c - M F )/{2M Bc ) ~ 0.8 GeV for decays to B s and ~ 0.9 GeV for de- 
cays to B mesons) . This means that the recoil momentum A of a final meson is of the same 
order as the relative momentum p of quarks inside a heavy-light meson (~ 0.5 GeV) in the 
whole kinematical range. Taking also into account that the final B s and B mesons are weakly 
bound [6], 1 we can replace the light quark energies by the center of mass energies on mass 
shell e q (p) -> E q = (Mj. -ml+ m 2 q )/(2M F ). We used such a substitution in our analysis of 
heavy-light meson mass spectra [6] which allowed us to treat the light quark relativistically 
without an unjustified expansion in inverse powers of its mass. Making these replacements 
and expansions we see that it is possible to integrate the current matrix element (11) either 
with respect to d 3 p or d 3 q using the quasipotential equation (1). Performing integrations 



The sum of constituent quark masses m& + m q is very close to the ground state meson mass M. 
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and taking the sum of the contributions and r^ 2 ^ we get the expression for the current 
matrix element, which contains ordinary overlap integrals of meson wave functions and is 
valid in the whole kinematical range. Hence the matrix element can be easily calculated 
using numerical wave functions found in our analysis of the meson mass spectra [5, 10]. 



IV. B r DECAY FORM FACTORS 



The matrix elements of the weak current J w for B c decays to pseudoscalar mesons (P = 
B s , B) can be parametrized by two invariant form factors: 



(P(pF)\qrc\B c (p Bc )) = f + (q 2 



Pb c + Pf 



Ml 



M 2 P 



Q 



+ fo(q 



2,M 2 Bc 



Ml 



q 



q 



(16) 



where q = p Bc — Pf] Mb c is the B c meson mass and Mp is the pseudoscalar meson mass. 

The corresponding matrix elements for B c decays to vector mesons (V = B*,B*) are 
parametrized by four form factors 



(V(p F )\qYc\B(p Bc )) = m^/I^ ^pPf,, 
(V(p F )\<ft li l5c\B( PBc )) = 2M v A (q 2 )^ q» + (M Bc + M v )A l (q 2 ) (e*» - ^ ^ 



(17) 



-A 2 (q 2 ) 



e -q 



Mi:. + My 



M Ml - My 

pb c +pf o — q 



(18) 



where My and e M are the mass and polarization vector of the final vector meson. The 
following relations hold for the form factors at the maximum recoil point of the final meson 
fa 2 = 0) 

/+(0) = /o(0), 

Mb - ~ Mv -Moy 



M0) = ^±Ml M0) 



■2 My V 7 2 My 

In the limit of vanishing lepton mass, the form factors f and A do not contribute to 
the semileptonic decay rates. However, they contribute to nonleptonic decay rates in the 
factorization approximation. 

It is convenient to consider B c semileptonic and nonleptonic decays in the B c meson rest 
frame. Then it is important to take into account the boost of the final meson wave function 
from the rest reference frame to the moving one with the recoil momentum A, given by 
Eq. (14). Now we can apply the method for calculating decay matrix elements described 
in the previous section. As it is argued above, the leading contributions arising from the 
vertex function can be exactly expressed through the overlap integrals of the meson 
wave functions in the whole kinematical range. For the subleading contribution T^, the 
expansion in powers of the ratio of the relative quark momentum p to heavy quark masses 
m^ c should be performed. Taking into account that the recoil momentum of the final meson 
A is not large we replace the final light quark energies e q (p) by the center of mass energies 
on mass shell E q . Such replacement is well justified near the point of zero recoil of the final 
B s , B meson. The weak dependence of this subleading contribution on the recoil momentum 
and its numerical smallness due to its proportionality to the small meson binding energy 
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TABLE I: Form factors of weak B c decays (c — > s(d) transitions). 



Transition 


7+(9 2 ) 


fo(q 2 ) 


V(q 2 ) 


Ai{q 2 ) 


A 2 {q 2 ) 


A (q 2 ) 


B c —> B S (B*) 














Q 2 = Q 2 


0.99 


0.86 


6.25 


0.76 


2.62 


0.91 


q 2 = 


0.50 


0.50 


3.44 


0.49 


2.19 


0.35 


B c -► 5(5*) 














a 2 = a 2 


0.96 


0.80 


8.91 


0.72 


2.83 


1.06 


g 2 = 


0.39 


0.39 


3.94 


0.42 


2.89 


0.20 



permits its extrapolation to the whole kinematical range. As a result, we get the following 
expressions for the B c decay form factors: 

(a) B c ^> P transitions (P = B s , B) 

fW) = /iV) +*/fV) + (1 -e)f V + ^{q\ (19) 

/o(g 2 ) = tfV) + eff 2 \q 2 ) + (1 - e)f^ 2 \ q % (20) 

(b) B c -> V transition (V = B*, B*) 

V(q 2 ) = V (1 \q 2 ) +eV S(2 \q 2 ) + (1 - e)V V(2 \q 2 ), (21) 

Mq 2 ) = A^(q 2 )+eA^ 2 \q 2 ) + (1 - e)A V ^\q% (22) 

A,{q 2 ) = 4\q 2 )+eA S 2 (2 \q 2 ) + (1 - e)A^ 2 \q\ (23) 

Ao(q 2 ) = 4V) + eA S (2 \q 2 ) + (1 - e)A^ 2 \q% (24) 

where /+q, f+^ 2 \ A^ 12 , Aq'J^ \ and V s ' v ^ are given in Appendix. The superscripts 
"(1)" and "(2)" correspond to Figs. 1 and 2, S and V to the scalar and vector potentials of 
gg-interaction. The mixing parameter of scalar and vector confining potentials e is fixed to 
be —1 in our model. 

It is easy to check that in the heavy quark limit the decay matrix elements (16)— (18) 
with form factors (19)-(24) satisfy the heavy quark spin symmetry relations [16] obtained 
near the zero recoil point (A — > 0). 

For numerical calculations we use the quasipotential wave functions of the B c meson, B s 
and B mesons obtained in the mass spectrum calculations [5, 6]. Our model predicts the B c 
meson mass Mb c = 6.270 GeV [10], while for B^ and fiW meson masses we use experimental 
data [17]. The calculated values of form factors at zero (g 2 = g^ax) an d maximum (q 2 = 0) 
recoil of the final meson are listed in Table I. In Fig. 3 we plot leading and subleading 
yS(2)^ jV(2) con t r i Du tions to the form factor / + for the B c — > B s transition, as an example. 
We see that the leading contribution f+* is dominant in the whole kinematical range, as it 
was expected. The subleading contributions f^ 2 \ are small and weakly depend on q 2 . 
The behavior of corresponding contributions to other form factors is similar. This supports 
our conjecture that the formulae (A.1)-(A.18) can be applied for the calculation of the form 
factors of B c — ► B S (B)^ transitions in the whole kinematical range. 
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In Figs. 4-7 we plot the calculated q 2 dependence of the weak form factors of Cabibbo- 
Kobayashi-Maskawa (CKM) favored (B c — > B s , B c — > B*), as well as CKM suppressed 
(B c — > B, B c — > 5*) transitions in the whole kinematical range. 

In the following sections we use the obtained form factors for the calculation of the 
semileptonic and nonleptonic B c decay rates. 



10 

6 
5 
4 
3 
2 
1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 

q 2 GeV 2 

FIG. 5: Form factors of the B c — > B*ev decay. 
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FIG. 6: Form factors of the B c — > Bev decay. 

V. SEMILEPTONIC DECAYS 

The differential semileptonic decay rates can be expressed in terms of the form factors as 
follows. 

(a) B c Pev decays (P = B s , B) 
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FIG. 7: Form factors of the B c — > B*ev decay, 
(b) B c -> 1/ez/ decays (V = B*, B*) 

- H = ^ (l"+te 2 )l 2 + l^-to 2 )! 2 + l^ofe 2 )! 2 

where GV is the Fermi constant, V cq is the CKM matrix element (q = s,d), 



A = IAI 



\ 



(Mi c + My -q 

AM 2 B 



2\2 



~ Mly. 



The helicity amplitudes are given by 

2M Bc A 



H±(q 2 ) 



M Bc + M v 



V(q 2 ) T iMBc + Mv) ! 1 ' • 



2M Bc A 



-AM 2 ) 



H (q 2 



2M V V¥ 



(M Bc + M V ){M 2 -M 2 - q 2 )Ai(q 2 ) 



4M| c A 2 



A 2 (q 2 



(26) 



(27) 
(28) 



M Bc + My 

The decay rates to the longitudinally and transversely polarized vector mesons are defined 
by 

\H (q 2 )\ 2 , 



dT L _ G 2 F A\V cq \ 2 q 2 



dq 2 



96vr 3 Ml 



dTV _ dr+ dT„ G 2 F A\V cq \ 2 q 2 
dq 2 dq 2 dq 2 



96tt 3 M% 



(\H + (q 



2\|2 



l^-(^)| 2 



(29) 
(30) 



In Figs. 8-11 we plot the differential semileptonic decay rates dT/dq 2 for semileptonic 
decays B c — > B^ev and B c — > B^ev calculated in our model using Eqs. (25), (26) both 
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FIG. 9: Differential decay rates (1/| V cs \ 2 )dT / dq 2 of B c -> B* s ev decay (in 1(T 12 GeV" 1 ). The 
upper curve is evaluated without account of l/mb iC corrections. 
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FIG. 10: Differential decay rate (1/| V cd \ 2 )dT/dq 2 of B c -> Bev decay (in 1(T 12 GeV" 1 ). The lower 
curve is evaluated without account of l/rrn, fi corrections. 
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FIG. 11: Differential decay rates (1/| V cd \ 2 )dT/dq 2 of B c -» B*ev decay (in 10~ 12 GeV" 1 ). The 
upper curve is evaluated without account of l/mb iC corrections. 
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TABLE II: Semileptonic decay rates V of B c to Bf' and BM mesons (in 10~ 15 GeV). 



Decay 


our 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


B c — > i? s ef 


12 


29 


59 


14.3 


26.6 


11.1(12.9) 


15 


12.3 


11.75 


26.8 


# c -> B^ezv 


25 


37 


65 


50.4 


44.0 


33.5(37.0) 


34 


19.0 


32.56 


34.6 


5c -> 5ef 


0.6 


2.1 


4.9 


1.14 


2.30 


0.9(1.0) 






0.59 


1.90 


£ c -► B*e^ 


1.7 


2.3 


8.5 


3.53 


3.32 


2.8(3.2) 






2.44 


2.34 



TABLE III: Semileptonic decay rates Tl^t,+ - (in 10 15 GeV) and their ratios for B c decays to 
vector B* and B* mesons. 



Decay 


T L 


r T 


r L /r T 


r+ 


r_ 


r+/r_ 


B c B*ev 


10.5 


14.5 


0.74 


3.1 


11.4 


0.27 


B c -► B*eu 


0.57 


1.13 


0.50 


0.13 


1.0 


0.13 



with and without account of l/m b c corrections to the decay form factors (A.1)-(A.18). 2 
From these plots we see that relativistic effects related to heavy quarks increase the rates 
of semileptonic B c decays to the pseudoscalar B s and B mesons and decrease the rates of 
semileptonic decays to vector B* and B* mesons. 

We calculate the total rates of the semileptonic B c decays to the B^ and B^ mesons 
integrating the corresponding differential decay rates over q 2 . For calculations we use the 
following values of the CKM matrix elements: \V CS \ = 0.974, \V c d\ = 0.223. The results are 
given in Table II in comparison with predictions of other approaches based on quark models 
[18, 20, 21, 23, 24, 26], QCD sum rules [19] and on the application of heavy quark symmetry 
relations [22, 25] to the quark model. Our predictions for the CKM favored semileptonic B c 
decays to B^ are smaller than those of QCD sum rules [19] and quark models [18, 20, 21, 26], 
but agree with quark model results [22-25]. For the CKM suppressed semileptonic decays 
of B c to B^ mesons our results are in agreement with the ones based on the application of 
heavy quark symmetry relations [22, 25] to the quark model. 

In Table III we present for completeness our predictions for the rates of the semileptonic 
B c decays to vector (B* and B*) mesons with longitudinal (L) or transverse (T) polarization 
and to the states with helicities A = ±1, as well as their ratios. 



VI. NONLEPTONIC DECAYS 



In the standard model nonleptonic B c decays are described by the effective Hamiltonian, 
obtained by integrating out the heavy W^-boson and top quark. For the case of c — > s,d 
transitions, one gets 



H cS = ^=V CS [cMO? + c 2 (i,)O c 2 s ] + -j=Vcd [cMOf + cMOf 



+ 



(31) 



2 Relativistic wave functions were used for both calculations. 
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The Wilson coefficients Ci j2 (/i) are evaluated perturbatively at the W scale and then are 
evolved down to the renormalization scale /i ~ m c by the renormalizat ion-group equations. 
The ellipsis denote the penguin operators, the Wilson coefficients of which are numerically 
much smaller than ci )2 . The local four-quark operators 0\ and 2 are given by 

Of = (du) V - A (cq)v-A, 

O2 = (cu) V - A (dq) V -A, q = (s,d), (32) 
where the rotated antiquark field is 

d = V ud d+V us s (33) 
and for the hadronic current the following notation is used 

(qq')v-A = ^(1 - l5)q' = J™ ■ 

The factorization approach, which is extensively used for the calculation of two-body 
nonleptonic decays, such as B c — > FM, assumes that the nonleptonic decay amplitude 
reduces to the product of a form factor and a decay constant [27]. This assumption in 
general cannot be exact. However, it is expected that factorization can hold for the energetic 
decays, where the final F meson is heavy and the M meson is light [28]. A justification of 
this assumption is usually based on the issue of color transparency [29] . In these decays the 
final hadrons are produced in the form of point-like color-singlet objects with a large relative 
momentum. And thus the hadronization of the decay products occurs after they are too 
far away for strongly interacting with each other. That provides the possibility to avoid the 
final state interaction. A more general treatment of factorization is given in Refs. [30, 31]. 

In this paper we consider the following two types of nonleptonic decays: (a) — > 
£«( J B(*)°)M+ and (b) B+ -> B^ + M°, where the final light M+ and M° mesons are tt, p 
or K(*\ The corresponding diagrams are shown in Fig. 12, where q, q\ = d, s and q 2 = u. 
Then in the factorization approximation the decay amplitudes can be expressed through the 
product of one-particle matrix elements 

(F°M + \H eS \B+) = ^V cq V qiq2ai (F\(cq) v ^ A \B c )(M\(q iq2 ) v _ A \0), 
(B^ + M°\H cS \B+) = ^V^V qiq MB { *\m2)v-A\B c )(M\(q iq ) v _ A \0), (34) 

where ^ ^ 

^ = ci(/i) + — c 2 (/i), a 2 = c 2 (/i) + — ci(/i) (35) 

and N c is the number of colors. 

The matrix element of the current between vacuum and final pseudoscalar (P) or 
vector {V) meson is parametrized by the decay constants fpy 

(^I?i7^75?2|0) = ifpft, (V%7^ 2 |0> = e,M v fy. (36) 

We use the following values of the decay constants: /„- = 0.131 GeV, f p = 0.208 GeV, 
f K = 0.160 GeV and f K * = 0.214 GeV. The CKM matrix elements are \V ud \ = 0.975, 
|K,|= 0.222. 
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b 



(a) 



(b) 



FIG. 12: Quark diagrams for the nonleptonic B c decays: (a) B+ -> F°M+ decay; (b) 5+ -> B+M° 
decay. 

The matrix elements of the weak current between the _B r meson and the final fiW, 5 
meson entering the factorized nonleptonic decay amplitude (34) are parametrized by the set 
of decay form factors defined in Eqs. (16) and (17). Using the form factor values obtained 
in Sec. IV, we get predictions for the nonleptonic £> + — > F°M + and — > B + M° decay 
rates and give them in Table IV in comparison with other calculations [19-23, 26]. 

In Tables II, IV we confront the predictions of our model for semileptonic and nonleptonic 
B c decays with previous calculations [18-26]. The constituent quark models of Refs. [18, 24] 
are based on the same effective quark-meson Lagrangian but use different phenomenological 
parameterizations (Gaussian [18] and dipole [24]) for the vertex functions, which are assumed 
to depend only on the loop momentum flowing through the vertex. The relativistic quark 
models of Refs. [20, 21, 26] use different reductions of the Bethe-Salpeter (BS) equation. 
The authors of Ref. [21] apply a nonrelativistic instantaneous approximation to the decay 
matrix elements and relate the BS wave functions to the Schrodinger ones. They do not 
give sufficient information about the quark interaction potential in their model. The quark 
models [20, 26] are based on the instantaneous approximation and different versions of the 
quasipotential equation. The meson wave functions are obtained by solving these equations 
with the one-gluon exchange plus the long-range scalar linear potentials. The light-front 
relativistic quark model is used in Ref. [23]. The decay form factors are expressed through 
the overlap integrals of the meson light-front wave functions which are related to the equal- 
time wave functions. The latter are expressed through the Gaussian functions. The heavy 
quark spin symmetry relations [16] and constituent quark models are used in Refs. [22, 25]. 
This symmetry permits to relate the B c weak decay form factors to a few invariant functions 
near the zero recoil point. These invariant functions are determined from the wave equation 
with the Richardson potential [22] or by the Gaussian wave functions [25]. Then they 
are extrapolated to the whole kinematical range accessible in B c decays. The authors of 
Ref. [19] employ three-point QCD sum rules with the account for the Coulomb-like a s /v- 
corrections. The values of the form factors are determined in the vicinity of q 2 = and then 
are extrapolated to the allowed kinematical region using the pole ansatz. 

Our relativistic quark model provides the selfconsistent dynamical approach for the calcu- 
lation of various meson properties. The meson wave functions in this approach are obtained 
as the solutions of the relativistic quasipotential equation. The weak decay matrix elements 
are expressed in terms of these wave functions. It allows us to determine explicitly the q 2 
dependence of the form factors of the weak B c decays in the whole kinematical range. All 
relevant relativistic effects are taken into account (including the boost of the wave functions 
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TABLE IV: Nonleptonic decay rates T (in 10" 15 GeV). 



Decay 


our 


[19] 


[20] 


[21] 


[23] 


[22] 


[26] 


B c 


— >■ i* s 7T 


25a? 


167a? 


-i f o 9 

15.8a? 


m a 9 

58.4a? 


34.8a? 


30.6a? 


^ K -I 9 

65.1a? 


B C 


D + 




72.5a? 


39.2a? 


44.8a? 


23.6a? 


13.6a? 


42.7a? 


B J 




16a? 


66.3a? 


-t r. r~ 9 

12.5a? 


51.6a? 


19.8a? 


35.6a? 


25.3a? 


B C 


o* + 


110a? 


204a? 


171a? 


150a? 


123a? 


110.1a? 


139.6a? 


B c 


— > t> s & 


2.1a? 


10.7a? 


1.70a? 


4.20a? 




2.15a? 


4.69a? 


B c ' 




0.03a? 




1.06a? 






0.043a? 


0.296a? 


B c 


— > B S K^ 


-i -I 9 

1.1a? 


3.8a? 


1.34a? 


2.96a? 




1.6a? 


1.34a? 


B t 


-^B°7T+ 


1.0a? 


10.6a? 


1.03a? 


3.30a? 


1.50a? 


1.97a? 


3.64a? 


B t 


-^B°p+ 


1.3a? 


9.7a? 


2.81a? 


5.97a? 


1.93a? 


„ _ . 9 

1.54a? 


. „ ^ 9 

4.03a? 


B t 




0.26a? 


9.5a? 


0.77a? 


2.90a? 


0.78a? 


. 9 

2.4a? 


1 nr. 9 

1.22a? 


B t 


5* V 


6.8a; 


26.1a? 


9.01a? 


11.9a? 


^ r-rr\ 9 

6.78a? 


8.6a? 


8.16a? 


B t 




0.09a? 


0.70a? 


0.105a? 


0.255a? 




0.14a? 


0.272a? 


B t- 




0.04a? 


0.15a? 


0.125a? 


0.180a? 




0.032a? 


0.0965a? 


B t- 




0.04a? 


n rr 2 

0.56af 


0.064a? 


0.195a? 




0.12a? 


0.0742af 


B+ -> B*°K*+ 


0.33a? 


0.59a? 


0.665a? 


0.374a? 




0.34a? 


0.378a? 


B t 




34a! 


286a! 


39.1a! 


96.5a! 


24.0a! 




103.4a! 


B t- 




13^ 


64a2 


/!/-> O 2 

46. 803 


68.2a! 


1 o o 2 

13.8a2 




36.6a! 


B t- 




15a! 


231a! 


24.0a! 


73.3a! 


8.9a! 




28.9a! 


B+ -> B*+K*° 


120a! 


242a! 


247a! 


141a! 


82.3a! 




143.6a! 


B t 




0.5a! 


5.3a! 


0.51a! 


1.65a! 


1.03a! 






B t 




0.7a! 


4.4a! 


1.40a! 


2.98a! 


1.28a! 






B t 


B*+TT° 


0.13a! 


4.8a! 


0.38a! 


1.45a! 


0.53a! 






B t 


-► 5*+/ 


3.4a! 


13.1a! 


4.50a! 


5.96a! 


4.56a! 







to the moving reference frame). The light quarks in the final heavy-light meson are treated 
relativistically. All that increases the reliability of the obtained results. 

As one sees from Tables II, IV the theoretical predictions for B c weak decay rates differ 
substantially. Thus experimental measurements of corresponding decay rates can discrimi- 
nate between various approaches. 

VII. CONCLUSIONS 

In this paper we considered weak semileptonic and nonleptonic B c decays to B s and 
B mesons, associated with c — > s,d quark transition, in the framework of the relativistic 
quark model based on the quasipotential approach in quantum field theory. The weak decay 
form factors were calculated explicitly in the whole kinematical range using the heavy quark 
expansion for the initial active quark c and spectator quark b. The final quark s or d 
was treated completely relativistically without applying an unjustified expansion in inverse 
powers of its mass. The leading order contribution of the heavy quark expansion was treated 
exactly, while in calculating the subleading order contribution the replacement of light quark 
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TABLE V: Branching fractions (in %) of exclusive B c decays calculated for the fixed values of the 
B c lifetime t Bc = 0.46 ps and a\ = 1.20, a 2 = -0.317. 



Decay 


or 


Decay 


r>r 


Decay 


or 


B c — ► B s ev 


U.84 


T~>-\- I ) 7>"" * + 


U.UU3 


T>-\- 7D*0 Ty*-\- 

BJ — > i? A T 


n no o 

U.U33 


±> c — > a s ev 


1.(0 


D+ , D* 7^ + 

£) c — ► fl s A 


nil 


R+ _^ R+ ^0 

_fc> c — > ±> A 




B c -> 


0.042 


B+ 5°7r+ 


0.10 




0.09 


5 C -► 5*61/ 


0.12 


S+ 5 V 


0.13 




0.11 


B+ _> B S 7T+ 


2.52 


5 + -► 5*°7T + 


0.026 


£+ B*+K*° 


0.84 


B+ - 


1.41 




0.68 


B + 5+vr 


0.004 


Bt - B*7T+ 


1.61 


B+ 


0.009 




0.005 


B+ - B>+ 


11.1 


B+ -► 5°^*+ 


0.004 


£+ B*+7r° 


0.001 




0.21 


5+ -► 5*°^+ 


0.004 


£+ B*+p° 


0.024 



energies e q (p) (q = s,d) by the center of mass energies E q on mass shell was performed. It 
was shown that such substitution introduces only minor errors which are of the same order 
as the higher order terms in the heavy quark expansion. Thus the decay form factors were 
evaluated up to the subleading order of the heavy quark expansion. The overall subleading 
contributions are small and weakly depend on the momentum transfer q 2 . 

We calculated semileptonic and nonleptonic (in factorization approximation) B c decay 
rates. Our predictions for the branching fractions are summarized in Table V, where we use 
the central experimental value of the B c meson lifetime [17]. From this table we see that the 
considered semileptonic decays to B s and B mesons give in total ~ 2.0% of the B c decay 
rate, while the energetic nonleptonic decays provide the dominant contribution ~ 19.3%. In 
our recent paper [3] we calculated weak B c decays to charmonium and D mesons, associated 
with b — > c,u quark transition. It was found that the semileptonic decays to the ground and 
first radially excited states of charmonium and to D mesons yield ~ 1.7% and corresponding 
energetic nonleptonic decays (to charmonium and K^> or ir,p mesons) contribute ~ 0.6%. 
All these decays (to B s , B, charmonium and D mesons) add up to ~ 23.6% of the B c total 
decay rate. 
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APPENDIX: FORM FACTORS OF WEAK B c DECAYS 

(a) B C ^P transition (P = B 8 , B) 



e q (p + A) + m q 



2e,(p + A) \ 2e c (p) 



e c (p) + m c 
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and the subscript g corresponds to s or d quark for the final B s or B meson, respectively. 



(b) B c -> V transition (V = 5;, 5*) 
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